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Abstract 

We investigate toy cosmological models in (1 + m + p)-dimensions with gas of p-branes wrapping 
over p-compact dimensions. In addition to winding modes, we consider the effects of momentum 
modes corresponding to small vibrations of branes and find that the extra dimensions are dynam- 
ically stabilized while the others expand. Adding matter, the compact volume may grow slowly 
depending on the equation of state. We also obtain solutions with winding and momentum modes 
where the observed space undergoes accelerated expansion. 
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I. INTRODUCTION 



String/M theory predicts six/seven extra dimensions which are presumably compact and 
very small compared to the observed three dimensions. To find out how extra dimensions are 
compactified is a major problem in the theory waiting for a solution. This is important in 
determining the vacuum structure and the low energy content which may then be compared 
with experimental findings to test validity of string/M theory. It is known that phenomeno- 
logically viable models can be obtained in Calabi-Yau or G2 manifold compactifications. 

Cosmology, on the other hand, is another framework where ideas in string/M theory may 
find observational support. We know that the perceived universe was initially small and 
grew to its size today after various cosmological eras. If extra dimensions exist they should 
have an effect on the cosmological evolution. However, since standard cosmology offers a 
successful scenario which remains plausible close to the temperatures around Tev scale, the 
effects of extra dimensions should be negligible from that time till today. For example, 
the theory of big-bang nucleosynthesis together with the observational abundances of light 
elements provide a limit to the change in the size of internal space following the production 
of these elements. Therefore, it is natural to assume that the extra dimensions were already 
stabilized around the epoch of nucleosynthesis. 

A mechanism for stabilization based on the cosmological impact of string winding and 
momentum modes was proposed in [1] (see also [2] and [3] for earlier work on string winding 
modes in cosmology). In the same paper, it is suggested that the initial big-bang singularity 
can be resolved by T-duality invariance of string theory and the dimensionality of observed 
space-time can be explained via the annihilation of winding strings leading decompactifica- 
tion of three dimensions. The proposal of [1] was generalized in [4] to include other branes 
in string/M theory and it is now an alternative approach named as brane gas cosmology (for 
recent work see e.g. [5-22]). In this scenario, strings dominate the the late time evolution 
and thus the stabilization of extra dimensions today should be achieved by string winding 
and momentum modes. This is quantified in the papers [19-22] where solutions to Einstein 
and dilaton gravity equations with stabilized internal dimensions are presented. 

However, this mechanism seems to work for toroidal compactifications which are not 
phenomenologically viable. For example, it is not clear how strings may stabilize volume 
of a Calabi-Yau or a G2 manifold since the topology is not equal to the product of one- 
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dimensional cycles (this is much more evident for a sphere compactification since there 
cannot be stable winding strings at all). Therefore, to fix the moduli related to higher 
dimensional non-trivial (i.e. non-toroidal) topological cycles, it seems that the higher di- 
mensional branes are also needed to play a cosmological role at late times. 

Motivated by this argument, in this paper we consider a toy cosmological model in (1 + 
m + p) -dimensions with a gas of p-branes wrapping over the p-dimensional compact space, 
and try to determine whether the compact volume is stabilized. We should emphasize that 
such a model does not straightforwardly lead to a realistic scenario. Apart from the fact that 
all other possible (stringy) excitations are ignored, the existence of a gas of 6 or 7-branes 
is also dubious (note that in a realistic model one should set m = 3 and p — 6, 7) since in 
the very early universe the higher dimensional branes are expected to annihilate leaving a 
gas of lower dimensional branes at late times. In any case, the results obtained in such toy 
models are still expected to be useful in brane gas cosmology. 

As in the original proposal of [1], one would anticipate a balance between brane winding 
and momentum pressures for stabilization. The energy momentum tensor for winding part 
was already determined in the cosmological context (see e.g. [9] or [14]). Since we do not 
know how to quantize p-branes for p > 1, the exact spectrum of the momentum modes is 
not known. However, in a semiclassical approximation one may consider small fluctuations 
around a rigid wrapped p-brane and from this spectrum an energy momentum tensor for 
momentum modes can be determined. In section II, we will work out this problem and 
study the issue of stability in the context of Einstein gravity. We find that, like strings, gas 
of p-branes with p > 1 can also dynamically stabilize extra dimensions. 

Based on recent observations one may claim that any potentially realistic cosmological 
model should include a theory for accelerated expansion. It is well known that inflationary 
paradigm is successful in resolving the basic shortcomings of standard cosmology. On the 
other hand, there are strong observational indications that the universe undergoes an ac- 
celerated expansion today. In brane gas cosmology, a way of realizing power-law inflation 
driven by a gas of domain walls was proposed in [23]. In section III, we also obtain ac- 
celerating solutions with brane winding and momentum modes (the background of winding 
modes generalizes the power-law metric constructed in [13]). As we will discuss, depending 
on the initial conditions, one encounters two different behavior: either a short period of 
acceleration with a small number (like 0.7) of e- foldings or an infinite accelerated expansion 
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ending at a singularity at finite proper time. 

II. STABILIZATION 

Consider a (1 + m + p)-dimensional space-time with the metric 

ds 2 = -dt 2 + dx i dx i + R 2 rfSj, (1) 

where i,j = l,..,m, S p is a p-dimensional compact manifold, R is a constant scale factor 
and the metric on S p can be written as (a, b — 1, ..,p), 

dY? p = g ab (y) dy a dy b . (2) 

In this paper we take g ab to be Ricci flat. For other cases, we expect the internal curvature 
terms to become important at late times as in [15]. In the physical gauge (and for small 
curvatures) the action for a p-brane wrapping over S p can be written as 




= -T p j d p+ V^ - ^ J cF+^y/^-yfdcx'dpx* + 0(x A ), (3) 

where T p is the brane tension, a a = (t, y a ) are the local brane coordinates, 7^ is the zeroth 
order induced metric given by 

-ylpdo" = -dt 2 + R 2 dT 2 pJ (4) 

and in the last line in (3) we expand the action to the second order in the transverse 
fluctuation fields x l (a a ). 

The first term in (3) corresponds to the winding modes which have the energy spectrum 

E = (nT p n p )R p , (5) 

where n is winding number and Q p is the volume of S p . From (5), the total pressure in the 
internal space can be found as 

P = -p(nT p Q p )RP, (6) 
which turns out to be negative, as expected. 
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On the other hand, (3) yields the following field equation for xf 

x i -±V 2 x i = 0, (7) 

where the dot denotes derivative with respect to t and V 2 is the Laplacian on £ p . Expanding 
in terms of the harmonics on E p , the spectrum of momentum modes (i.e. the eigenfrequencies 
e luJt ) can be obtained as uj 2 = X^/R 2 , where —A 2 are the eigenvalues of the Laplacian. Thus 
the energy for the n'th momentum mode is given by (we choose A n > 0) 

E = ^ (8) 
and the corresponding total pressure on S p is 

P=% (9) 

which is positive contrary to (6). Although (8) and (9) depend on the undetermined eigen- 
values A n , one does not need these numbers to proceed. On the other hand, 1/R dependence 
in (9) can also be deduced from the uncertainty principle Ay AP ~ 1 where a compact di- 
rection y in S p approximately satisfies Ay ~ R. This suggests that (8) and (9) should also 
hold for other possible fields propagating on the world- volume such as U(l) gauge fields on 
D-branes or self-dual 3-form of M5 brane. Gauging away the longitudinal polarizations, the 
linearized transverse modes of these brane fields are expected to obey (7). 
In a cosmological setting, (1) should be replaced with 

ds 2 = -e 2A dt 2 + e 2B dx l dx l + e 2C dT 2 , (10) 

where the metric functions A, B and C depend only on time t. The scale factors for the 
observed and the internal spaces are defined as 

R ob = e B , R m = e c . (11) 

Note that t is not the proper time and we have not yet fixed the corresponding reparametriza- 
tion invar iance. 

From (5)- (9), the energy momentum tensor for the brane winding and momentum modes 
can be obtained by calculating the corresponding densities (i.e. by dividing with the total 
spatial volume V = e mB+p ) which gives 



rp _ rp —mB _|_ rp —mB—(p+l)C 
± tt — ± w e t -t m t- j 



T Vj = 0, (12) 
T hl = -T w e- mB 5 ab + ^e- mB ^ c 5 ab , 

where T w and T m are constants (which can in principle be fixed by the parameters of the 
underlying fundamental theory and thermodynamics) and the hatted indices refer to the 
tangent space. It is easy to verify that V^T^ = 0. In driving (12), the total pressures given 
in (6) and (9) are assumed to be distributed isotropically on the tangent space of E p since 
the momentum modes are point-like excitations (note that we are doing field theory on the 
world- volume) and thus locally they would only see the flat, tangent space. However, to 
justify this assumption more rigorously one may additionally require that S p is not highly 
curved and thus has a "smooth" shape with no sharp corners or handles, which is also 
necessary for the validity of the brane action and the small field approximation used in (3). 
Let us note that the contribution of the winding part in (12) agrees with the expressions 
derived in [13, 14] which are obtained in a different way by coupling the brane action to the 
gravity action. This agreement further supports the above considerations. 
In general, if one has the metric 

ds 2 = -e 2A dt 2 + ]T e^dx'dx', (13) 

and the energy momentum tensor 

T A » = diag( M ), (14) 
where p- = uj{ p with constant u>i, the conservation equation V M T^ = gives 



p = po exp 



(15) 



where p is a constant. Comparing to (12), we find that the energy momentum tensor for 
winding and momentum modes can be written as (14) where for winding modes 



— 1 : brane direction, 

(16) 

: transverse direction, 



and for momentum modes 

1/p : brane direction, 



(17) 

: transverse direction, 
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which corresponds to radiation confined in the compact space. 1 With the above formulas it 
is possible to determine the energy momentum tensor for any given brane configuration and 
study the resulting cosmological evolution. 

Using (12) and imposing the gauge A = mB + pC in (10), which fixes the t- 
reparametrization invariance, Einstein equations can be written as follows: 

A - A 2 + m B 2 + p C 2 = m ~ 2 T w e mB+2pC - T m e mB+ ^ c , 

m + p — 1 

B = P+1 T w e mB+2pC , (18) 
m + p — 1 

r — - 771 ~ ^ rp mB+2pC _|_ mB+(p-l)C 

m + p — 1 p 
where the dot denotes derivative with respect to t and the gravitational constant is set to 
one. 

We could not obtain the most general solution of (18). However, it is easy to verify that 
the functions 

_2 

A = -21n(at) +pC , B = — \n(at), C = C , (19) 

m 

give a solution where a is a constant, which is uniquely 2 fixed by field equations, and Co 
can be calculated from the third equation as 



in 



, T w (m-2)p _ 



(20) 



Introducing the proper time coordinate 

dr = -e A dt, (21) 

we get 

ds 2 = -dr 2 + [arfl m dx l dx l + (R° in ) 2 <E 2 , (22) 

where the internal dimensions are stabilized and the power-law of the observed space is 
exactly the same as the one for pressureless matter in standard cosmology. Here, R® n can 



1 The energy-momentum tensor for (16) or (17) can also be viewed to correspond to a perfect fluid in 
(1 + m + p)-dimensions. Certain aspects of higher dimensional cosmological models with general multi- 
component perfect fluid matter were studied in the past, see e.g. [24-27]. 

2 The field equations give 

2 m{p+l)T rn j 
a ~ 2p(m-2) {Rin) • 
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be viewed as the self-dual radius corresponding to brane winding and momentum modes. 
Technically, C is allowed to be a constant since in the third equation in (18), exp(mB) terms 
can be factored out and there are both negative and positive contributions to C. 

The metric (22) corresponds to the initial data where the internal dimensions begin at 
the self-dual radius with zero velocity. For more general initial conditions one may integrate 
(18) numerically. The first equation (when combined with the other two) in (18) gives a 
constraint on the initial data: 

(23) 

where r is the proper time defined in (21). Fixing B(t), C(t) and one of the velocities 3 
dB/dr or dC/dr at some r = r , (23) gives two real roots for the undetermined initial 
velocity. It is easy to see that one of the roots should have the opposite sign compared to 
the other fixed initial velocity. 

After several numerical integrations of (18) for different initial conditions and for 
(m,p) = (3,6) (these (m,p) values are suggested by string theory, however we expect the 
same conclusions to hold for any set with m > 2) we observe the following generic behavior: 

(i) For one of the roots, exp(C) = Ri n performs damped oscillations around the self-dual 
radius given in (20) and the observed space expands with negative acceleration close to 
the power-law in (22). As the extra dimensions are stabilized asymptotically, the solution 
becomes (22). In figure 1, we plot two such numerical runs for R in . As in [19], the evolution 
of the large dimensions is responsible for these oscillations to be damped. 

(ii) For the other root, the numerical integration signals a singularity at finite proper time, 
where the observed and the internal spaces have the opposite behavior near the singularity 
i.e. if one expands the other contracts. Below, we will try to explain this case analytically. 

(iii) If the initial values of dB / dr and dC / dr are positive (i.e. if both spaces are expanding 
at the beginning) then we observe the stabilization mentioned in (i). For the other cases 
both (i) and (ii) are possible. 

It is not surprising that for some initial conditions the metric ends with a singularity at 
finite proper time. We do not expect to obtain a completely smooth cosmological solution 
with the ansatz (10). Contrary, at least one singular fixed point in the past or in the future 

3 In our numerical integrations, we first fix dC/dr and solve dB/dr from (23). 
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FIG. 1: The graphs of Ri n {r) for two different initial data set obeying (i). For convenience we 
take m = 3, p = 6, T w = 8 and T m = 6 so that the self-dual radius (20) is equal to 1. The initial 
conditions are: R ob (0) = 1, R in (0) = 0.7, dR ob /dT(0) = 3.9761, dR in /dT{0) = 0.3 and R ob (0) = 1, 
R in (0) = 1.3, dR ob /dr(0) = 3.65135, dR in /dT(0) = -0.7. Here, dR ob /dT(0) is solved from (23) 
given other conditions. 

is anticipated. Therefore, for the alternatives (i) and (ii) the following interpretation can 
be given: in the phase space one of the roots selects a trajectory which evolves "forwards 
in time" and gives the smooth stabilized solution, and the other root corresponds to a path 
where the space-time evolves "backwards in time" and hits the initial singularity. 

Actually, to define the arrow of time in such solutions an extra (observational) input 
is needed. For instance, we use the Robert son- Walker metric starting from the big-bang 
singularity since we observe that the universe at the moment is expanding. If we would 
measure a contraction rather than an expansion, the same metric should be considered with 
r — > — t, i.e. with a future big-crunch singularity instead. In our case, there is no such 
avaliable observational input. However, it seems reasonable to assume that (18) holds some 
time after the big-bang explosion and thus one should take all dimensions to be expanding 
initially in (18). From (iii) we then see that the corresponding solution asymptotes to (22) 
and the stabilization can be achieved. On the other hand, we should also bear in mind that 
there are physically acceptable initial data, such as an expanding observed and contracting 
internal dimensions at the beginning, where the subsequent evolution ends with a naked 
singularity at finite proper time. 

To have a better understanding of this singular behavior, let us analyze (18) more care- 
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fully. Comparing the contributions of the winding and momentum modes, we see that when 
exp[(p+l)C] 3> T m /T w the winding modes and when exp[(p+l)C] <C T m /T w the momentum 
modes dominate the evolution. So, it would be useful to solve (18) in these two limits. 

To consider the effects of winding modes alone we set T m = 0. Then the last two equations 
in (18) implies 

m — 2 

C = -—-^B+at, (24) 
p+1 V ' 

where a is a constant. Using this relation in (18), one gets a second order differential equation 
for B which can be solved (up to an irrelevant additive constant depending on T w ) as 

B = - 2(P+1 \ ln(sinh(6t)) - 2P{P+1 \ at (25) 
m — mp + 4p m — mp + 4p 

where b is another constant. Using (24), (25) and the gauge A = m B + pC in the first 

equation in (18) one finds 

B , = 4(m + p-l ) t , 
mp(p + iy 

and this completes the solution. In terms of the proper time coordinate (21), we see that as 
t — > 0, r — > oo and as t — > oo, r — > 0. In these two limits the solution becomes 

r -dr 2 + (r) 2ri rf^ rfx* + (r) 2r2 rfS 2 : as r -> 0, 

^4 = (27) 

4 4(m-2) 

[ -rfr 2 + (r)™ rfx 1 rfx 1 + (r) rfS 2 : as r -> oo, 

where the powers ri and r 2 are given by 



m + e \/mp (m + p — 1) 

ri = ( M , -» (28) 

m(m + p) 



p — e \/mp(m + p — 1) 

^ = V , ^ -■ (29) 

p(m + p) 

Here, e = ±1 corresponds to the positive and the negative roots of (26), respectively, and in 
(27) we ignore constants multiplying r in parenthesis which can be set to 1 by scalings of 
x and y coordinates. We see that the full metric smoothly interpolates between these two 
limiting geometries. 

Asymptotically, as r — > oo the solution approaches the power-law metric given in [13] 
(see eq. (11) in that paper). On the other hand, as r — > one finds a vacuum Kasner metric 
since the powers r\ and r 2 obey 

mr\ + pr 2 = 1, (30) 
mr\ + pr\ = 1. (31) 
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Indeed, using (24) and (25) it is easy to see that all the source terms on the right hand 
side of (18) vanish as t — > oo (and thus r — > 0) so ending up with a vacuum solution is not 
surprising. This also shows that the same vacuum metric solves (18) with T m ^ in the 
t — > limit. 

To see the effects of momentum modes alone let us now take T w = 0. In this case the 
second equation in (18) gives 

B = at, (32) 

where a is a constant, and the third equation can be solved for C (up to an irrelevant additive 
constant depending on T m ) as 

2 m 

C = ln(sinh(6t)) a t. (33) 

p — 1 p — 1 

The unknown function A can be fixed from the gauge condition A = mB + pC, and the 
first equation in (18) imposes 

a 2 = —f-^- ~rb 2 . (34) 

m{m + p — 1) 

In terms of the proper time (21), we find that as t — > (0, oo), r — > (oo, 0), respectively, and 
asymptotically the solution becomes 

' -dr 2 + (r) 2ri dx i dx i + (r) 2ra dY? p : as r -> 0, 
tfc& = (35) 

4 

where ri and r 2 are given by (28) and (29) in which e = +1 corresponds to the negative root 
of (34) while e = — 1 corresponds to the positive root. As r — > 0, one again encounters the 
same vacuum Kasner metric. 

Returning back the situation where both T w and T m are non-zero, we see from (27) and 
(35) that as r — > the solution should approach the Kasner background with the powers 
T\ and r 2 . On the other hand, the field equations show that as r — > oo the metrics given in 
(27) and (35) should be modified i.e. in (27) momentum modes can no longer be ignored 
and in (35) winding modes should be taken into account. In either case the neglected 
contribution becomes dominant in time, alters the behavior of the extra dimensions and 
serves for stabilization. This suggests that asymptotically the most general solution of (18) 
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should become 

( -dr 2 + (r) 2ri dx i dx { + (r) 2r2 dY? p : as r -> 0, 
rfs 2 = | (36) 
[ -dr 2 + (r)£ tfo* dx* + (R° in ) 2 dY? v : as r -> oo. 

The numerical integrations of (18), which are summarized in (i)-(iii) above, also support this 
result. We already mentioned in (i) that for some initial conditions the metric approaches 
to (22). The singularity corresponding to other initial data indicated in (ii) occurs at r = 
in (36). Recall from (ii) that in reaching the singularity the observed and the internal spaces 
have the opposite behavior. From (28) and (29), we see that r\ and r 2 have the opposite 
signs and thus (36) also explains this numerical result. As discussed above, such initial 
conditions are associated with solution evolving "backwards in time" . 

Having studied the stabilization problem with brane winding and momentum modes, let 
us now try to include the effects of ordinary matter. These should be taken into account 
in determining the late time behavior since in the universe today the existence of radiation 
and pressureless dust is known. We take the following energy momentum tensor for matter 

T Al > = diag(p, p h pa), (37) 

where 

p-=ujp, Pa = vp, (38) 

and (u,v) are constants. Adding (37), the Einstein equations (18) are modified as follows 
(we again impose the gauge A = m B + pC): 



2 - m(u + 1) -p{y + 1) 
m + p — 1 



•• ' ' Tft — 2 

A- A 2 + mB 2 +pC 2 = ■ T W F W — T m F m + p 

m + p — 1 

B= P A + 1 - T W F W + Po [l + (p-l)u-pu]F p , (39) 
m + p — 1 m + p — 1 

~ = m-2 TwFw + Trn Fm + Po [i + (m _i) z/ _ mCL; ] j p 

m + p — 1 p m + p—1 



where 



P . mB+2pC 



mB+{p-l)C 



P = e {l-w)mB+(l-v)pC _ 



(40) 
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It is difficult to solve (39) in its most general form. However, for uj = 0, there is again a 
special solution which can be written in the proper time coordinate as: 

ds 2 = -dr 2 + {aT) A/m dx l dx l + (i?° n ) 2 d£j, (41) 

where the self-dual radius i?° n obeys 



m 



~~~r T w (Rir 1 - 1 + (m l ] V po {Rlf- vp - T - = 0, (42) 



m + p — 1 m + p — 1 p 

and a is a fixed 4 , positive constant. 

The metric (41) evolves from the initial data where the internal space starts at the self- 
dual radius with zero velocity. For more general initial conditions, we numerically integrate 5 
(39) and our findings can be summarized as follows: 

(a) uj = 0: We find that the internal dimensions are stabilized after some damped os- 
cillations around the self-dual radius determined in (42) and the observed space expands. 
Asymptotically, the solution approaches to (41). In figure 2(a), two such numerical inte- 
grations for Rin(r) are plotted. Note that, from (42) the final size of the extra dimensions 
depends on v. 

(b) uj > 0: We find that the effect of ordinary matter on the evolution is negligible at 
late times, i.e. brane winding and momentum modes dominate the cosmology. The solution 
asymptotically becomes (22) where the self-dual radius given by (20) does not depend on v. 
In figure 2(b), we plot two numerical runs for R in {r). The initial conditions are chosen so 
that R° in = 1. 

(c) uj < 0: In this case, the late time evolution is dominated by matter and the stabiliza- 
tion cannot be achieved. Two numerical runs for Ri n {r) are plotted in figure 2(c). At early 
times, the oscillations caused by winding and momentum modes can be seen in the figure. 
However, matter finally causes the internal space to grow at late times. 



The field equations impose 



a 2 



m(p + l)T m x m{v + l)p a p 
2p(m-2) { in> + 2(m-2) ( in) 



5 As before, we first fix all but the velocity of the observed space at t = and solve dB(r)/dT from the 
constraint equation which gives two real roots. In the following, we only consider one of the roots which 
yields an evolution "forwards in time" . 
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(c) 

FIG. 2: The graphs of Ri n {r) for different values of to and v. For convenience we take m = 3, 
p = 6, T w = 8, T m = 6, R o b(0) = 2. The initial conditions are: 

(a) [uj = 0]: v = -0.9, i? in (0) = 1.1, dR ob / dr{<d) = -1.41365, dR in /dT(0) = 1; 
v = 0.5, R in {0) = 1.1, dR ob /dT(0) = 6.80001, dR in /dT(0) = -0.7, 

(b) [w = 0.5]: i/ = -0.6, R in {0) = 1.01, dR ob /dr{{)) = 1.10393, dR in /dT(0) = 0.1; 
i/ = 0.3, i? in (0) = 0.8, dR ob /dr(0) = 1.07248, dR in /dr(0) = 0.3, 

(c) [w = -0.1]: v = -0.3, J2i„(0) = 1.7, dR ob /dT(0) = -0.870317, dR in /dT(0) = 1; 
z> = 0.4, i? in (0) = 1.1, dR ob /dT(0) = -1.26295, dR in /dr{{)) = 0.9. 

Here, dR ob /dr(0) is solved from the constraint on initial data. 

From (39), we see that the contributions of matter, the winding and the momentum modes 
are proportional to the functions F p , F w and F m given in (40), respectively. Assuming that 
the size of the observed space becomes much larger than the size of the extra dimensions in 
time 6 , one sees that e B factors in (40) dictate a hierarchy between F- functions. This explains 

6 This is expected since the winding modes resist the expansion of the extra dimensions, however all modes 
contribute to the growth of the observed space. 
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the generic behavior we encountered in (a)-(c): compared to F w and F m , the function F p 
gets stronger for u < and it becomes negligible for uo > 0. In the former case it is clear 
that the stabilization mechanism does not work. When uo = all functions have the same 
strength and stabilization can be achieved due to the presence of positive and negative 
contributions to C in (39), where the self-dual radius now depends on the equation of state. 
These considerations suggest that if the acceleration of the universe is due to some form 
of matter with uo < (i.e. dark energy), then the extra dimensions cannot be stabilized 
by brane winding and momentum modes. We expect the same conclusion to hold for the 
strings studied in [19]. However, depending on the actual values of the parameters in these 
equations, the change in the size of the extra dimensions may be slow enough for us not to 
recognize it in our time scale. 



III. ACCELERATION 



In this section, we show that some of the solutions constructed above with T rn = or 
T w = have periods during which the observed space undergoes accelerated expansion. 
Here, we do not intend to propose a realistic model for inflation or current acceleration 
of the universe, but simply note an alternative way of obtaining acceleration in brane gas 
cosmology. Motivated by string/M theory as the natural fundamental framework, in this 
section we also set m = 3 and p = 6 for convenience. 

Let us consider the solution for the brane winding modes given in (24)- (25). Setting b = 1 
and choosing the negative root 7 in (26), the metric functions become 

10 / ■ i \ 8 
A = — — ln(smht) + - t, 
o o 

14 16 
B = -— ln(sinht) + — t, (43) 

C = -ln(sinht)-^, 

where the additive constants, which can be shifted by scaling out x and y coordinates in 
(10), are set to zero. Eq. (21) can be integrated to fix the proper time: 

r = — [cosh(t/3) - sinh(t/3)] [-21 + 23 cosh(2t) + 5 sinh(2t)] sinh(t)" 7/3 . (44) 



7 For the positive root, we find that there is no period of accelerated expansion. 
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FIG. 3: The graphs of v b in (45) and a D b in (46) with respect to the proper time r in (44). 

As pointed out before, in the limits t — > 0, oo we have r — > oo, 0, respectively. Using (43), 
it is easy to calculate the proper velocity and the acceleration of the observed space 

dR b 2 



Vob = 
Clob = 



dr 9 
oh 2 
81 



fR 

~dV 2 



e -8t/9 [_g + 7coth (f)] sin h(t) 16 /9 ? (45) 

e -32t/9 [_ 3g + 8g cosh (2t) _ 92 sinh(2t)] sinh(t) 28/9 , (46) 



where the scale factors R Q b and R in are defined in (11). In figure 3, we plot v Q b and a Q b 
versus proper time r. 

From (44), (45) and (46) one finds that v Q b = at r = T\ ~ 6.5 and a ob = at 
t = t 2 ~ 23.3. Moreover, i> o6 is positive for r > T\ and a ob is positive for r < r 2 (see figure 
3). Therefore, in the interval T\ < r < r 2 we have both v Q b > and a ob > which implies 
an accelerated expansion for the observed space. The corresponding number of e-foldings is 
approximately 0.7. 

It is clear that this growth is not enough for early inflation, however the metric may 
explain the current observational acceleration of the universe. In the solution, the internal 
space expands till r = r 3 ~ 9.9 and then it begins to contract. It is interesting to note that 
r 3 G (ti, r 2 ) and thus acceleration occurs around when the internal space alters its evolution 
from expansion to contraction. 

On the other hand, it is also possible to view the metric "backwards in time" by sending 
t — > —t. Although this does not change the acceleration, the velocity picks up a sign. In the 
reparametrized solution, we see that as r runs from —T\ to we have v ob > and a Q b > 0. 
This is an infinite accelerated expansion in a finite proper time interval and the metric ends 
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with a singularity at r = 0. As r — > , the solution approaches to 

ds 2 = -dt 2 + (-ry^dx'dx' + {-rf^dT 2 , (47) 

which is the time reversal of the vacuum Kasner metric obtained in (27) with e = — 1 and in 
(35) with e = l. Note that the acceleration observed in (47) is different than the power-law 
inflation and it would be interesting to calculate the spectrum of primordial fluctuations to 
compare the results with the measured spectral indices. 

Accelerating cosmologies can also be obtained with brane momentum modes. After send- 
ing t — > — t, the metric functions for the solution (32)- (33) become 

12 3 
A = ln(- sinh t) + -t, 

5 5 

B = -t, (48) 

2 3 
C = -- ln(-sinht) + -t, 

5 5 

where in (34) the positive root 8 is picked up and we again set b = 1. Note that t is defined 
in the negative line (— oo, 0). It is possible to choose the integration constant in (21) so that 
as t — > — oo, we have r — > — oo, 0, respectively, and thus r is also defined in the negative 
line. Calculating the proper velocity and the acceleration of the observed space one finds 

dR n 



Vob = 



dr 



^ = e- 8 '/ 5 sinh(-t) 12 / 5 , (49) 



_ d 2 R ob _ nt/5 r 12 
a ob - ^ 2 



cosh(— t) H — sinh(— t) 

5 5 



sinh(-t) 19 / 5 , (50) 



which are always positive and thus give an accelerated expansion. Correspondingly, one can 
see that the internal space monotonically contracts. This evolution ends with a singularity 
at t = where the metric becomes (47) as r — > 0~. The expansion is limited near the 
asymptotic region r — > — oo (i.e. as t — > — oo) since both t> ob and a ob vanish. However, a 
large number of e-foldings can be obtained by using the solution near r = 0. 

IV. CONCLUSIONS 

In this work, based on the cosmological impact of brane winding and momentum modes, 
a stabilization mechanism for extra dimensions is proposed. We find that a gas of p-branes 



It turns out that the negative root does not yield acceleration. 
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wrapping the p-dimensional compact space can stabilize its volume as a result of a balance 
between winding and momentum pressures at the "self-dual" radius. This generalizes the 
method studied in [19] to p-branes with p > 1. As discussed in the introduction, a mechanism 
involving p-branes is desirable since strings cannot keep the volume of a generic internal 
manifold constant for topological reasons. On the other hand, one should also bear in mind 
that it is difficult to build late time realistic scenarios involving p-brane gases with p > 2, 
since the higher dimensional branes are likely to annihilate in the early universe. 

Adding matter, we find that the fate of the stabilization depends on the equation of state 
p = up in the observed space. For uj > 0, the matter can be ignored at late times and for 
u) = the self-dual radius is modified but stabilization can still be achieved. When uj < 0, 
the late time dynamics is dominated by matter and one finds an expanding internal space, 
which shows that the stabilization mechanism does not work. 

In this paper, we also construct accelerating solutions driven by brane winding and mo- 
mentum modes. For winding modes, a short period of acceleration can be observed when 
the evolution of the internal space changes from expansion to contraction. This suggests an 
interesting mechanism for acceleration which deserves further study. Since the number of 
e-foldings for the expansion is small, this metric cannot be used to explain the early inflation. 

In some solutions there are also accelerating periods ending with a naked singularity, 
which can be seen both for winding and momentum modes. Here, a large number of e- 
foldings can be obtained as one approaches the singularity, where the solution becomes one 
of the vacuum Kasner metrics. Such backgrounds may be useful in realizing inflation in 
brane gas cosmology, however it is necessary to find a way of resolving the singularity first, 
which may also explain the exit from inflation. Since the internal space shrinks to zero size, 
the momentum modes are expected to become important. Recall that in some solutions 
the metric evolves through the singularity in the presence of linearized momentum modes. 
However, the small-field approximation should brake down near the singularity and the full 
momentum modes might stop the contraction. T-duality invariance of string theory may 
also play a crucial role here, so it would be interesting to include the dilaton into the system. 
In any case, the quantum gravitational effects are expected to cure the singularity, so one 
can try to cut the solution and paste it to some other metric to describe the subsequent 
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evolution. Studying intersecting branes might also be important in this context. 
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